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In this paper we deal with neutron stars, which are described by a perfect fluid model, in the context of the 
teleparallel equivalent of general relativity. We use numerical computations (based on the RNS code) to find 
the relationship between the angular momentum of the field and the angular momentum of the source. Such 
a relation was established for each stable star reached by the numerical computation once the code is fed with 
an equation of state, the central energy density and the ratio between polar and equatorial radii. We also find a 
regime where linear relation between gravitational angular momentum and moment of inertia (as well as angular 
velocity of the fluid) is valid. We give the spatial distribution of the gravitational energy and show that it has a 
linear dependence with the squared angular velocity of the source. 
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I. INTRODUCTION 

Most of the known stars are accurately described by New- 
tonian physics. Such compact objects as neutron stars with 
masses of the order of one solar mass, by contrast, are not. 
The gravitational field of such objects is so strong that only 
general relativity is able to explain their properties [1 ]. Given 
that, due to the very nature of Einstein's equations, analytical 
methods are very difficult to implement, most interesting re- 
sults come from non-algebraic methods. In particular, to deal 
with with rotating neutron stars, one usually resorts to numer- 
ical procedures ||2|,|3|]. 

While the existence of such exotic objects was theoretically 
predicted in the 1930's fllD], the first observations took place 
in the 1960's J^-IH]. In neutron stars, the mechanism pre- 
venting gravitational collapse is a repulsive interaction due to 
the quantum nature of the particles in the star, which com- 
pensates the gravitational pressure, so that neutron stars exist 
in very compact spatial volumes. Neutron stars with masses 
are around one solar mass and radii are around 10 Km are 
known [9], as are pulsars rotating with periods as low as a 
millisecond and very high surface magnetic fields ifToL [TTll . 
These very unusual conditions make them perfect candidates 
to test any new predictions of a gravitational theory. 

We intend to treat rapidly rotating neutron stars from the 
viewpoint of teleparallel gravity. From the dynamical point of 
view, teleparallel gravity and general relativity make identical 
predictions. On the other hand, teleparallel gravity allows for 
the definition of physically interesting quantities, such as the 
the well-behaved gravitational field energy-momentum and 
angular- momentum tensors 01211 . 

In general relativity, the definitions of these quantities are 
still under development. The first attempt to define gravita- 
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tional energy was based on pseudo-tensors that were not in- 
variant under coordinate transformations. This effort was fol- 
lowed by the Komar integrals [13] and then by the formalism 
known as the ADM formulation 1 3, which is based on the 
(3+l)-dimensional Hamiltonian formulation of general rela- 
tivity, an approach using constraints to define energy, momen- 
tum and angular momentum. The difficulty with this formal- 
ism is that such quantities are only well-defined asymptoti- 
cally, at the spatial infinity. In addition, no expression ob- 
tained in the realm of general relativity is dependent on the 
reference frame — clearly a feature that is undesirable for en- 
ergy, momentum, and angular momentum. 

The expressions for the energy-momentum and angular mo- 
mentum of the gravitational field, in the context of the Telepar- 
allelism Equivalent to General Relativity (TEGR), are invari- 
ant under transformations of the coordinates of the three- 
dimensional spacelike surface. They are also dependent on the 
frame of reference, as one would expect. Over the years, they 
have been consistently applied to many different systems IU5l - 

a. 

The frame dependence is an expected condition for any ex- 
pression due to the field since in special relativity the energy of 
a particle for a stationary observer is m, but it is ym for an ob- 
server under a Lorentz boost, where c = 1 and y is the Lorentz 
factor. There is no reason to abandon this feature when deal- 
ing with the gravitational field. Similar considerations apply 
to the momentum and angular momentum. 

Our goal is to establish a relation between the angular mo- 
mentum of the gravitational field, which is only predicted by 
the teleparallel theory of gravity, and the angular momentum, 
the angular velocity of the source, and the moment of iner- 
tia. We hope to identify the significance of the field quantities 
when related to matter quantities. In particular, it is known 
that in the final stage of the coalescence between two black 
holes in a binary system the remnant black hole will acquire 
linear momentum As a matter of fact, the field has to 

gain the exact amount needed to cancel it in order to preserve 
momentum conservation [22]. This example highlights the 



Ulhoa and Rocha 



2 



importance of knowing the relation between field and source. 
We also intend to describe the behavior of the energy over the 
entire spacetime and to relate it to other features of neutron 
stars. 

This paper is organized as follows. In Sec. [II] we briefly 
discuss certain ideas in teleparallel gravity. In Sec. |III]we de- 
scribe the space-time of a neutron star and fix the frame that 
will be used to compare features only due to the gravitational 
field with others due to the source. In Sec. lIII Al we discuss the 
field angular momentum and its relation to certain features of 
the source. In Sec. IIIIBl we relate the gravitational energy to 
the angular velocity of the fluid. Plots of the numerical results 
are presented to describe the distribution of gravitational en- 
ergy in space. Finally, the last section contains our concluding 
remarks. 

In our notation, space-time indices /j,v, ... and SO(3,l) in- 
dices a,b, ... run from to 3. Time and space indices are in- 
dicated as /j = 0,z", a = (0), (/). The metric tensor g pV raises 
and lowers space-time indices, while the Minkowski metric, 
\\ ab = diag{ - + ++), acts on SO(3,l) indices. The tetrad 
field is denoted by e a », and the determinant of the tetrad field 
is represented by e = det(e" /J ). Unless otherwise stated, we 
adopt units such that G = c = 1 . 



II. TELEPARALLELISM EQUIVALENT TO GENERAL 
RELATIVITY (TEGR) 

Teleparallel gravity has been investigated over the years as 
an alternative to the theory of general relativity (GR). Such a 
theory is entirely equivalent to GR as far as the field equa- 
tions are concerned. Both theories are derived from La- 
grangians that only differ by a total divergence, and therefore 
have the same field equations. The main advantage of TEGR 
over GR is allowing us to construct consistent expressions for 
energy, momentum, and angular momentum 01211 . We will 
now present ideas and expressions developed in the realm of 
TEGR. 

The dynamical variables of TEGR are the tetrad fields de- 
fined in the Weitzenbock space-time, also known as the Cartan 
space-time. In the Riemann space-time the dynamics is given 
by the metric tensor. The tetrad field and metric tensor are 
related by the equalities 

g^^eV. (1) 

Therefore, for each metric tensor one can define infinitely 
many tetrads. It follows that the Weitzenbock geometry is 
less restrictive than Riemann geometry. The arbitrariness in 
the choice of the tetrad field is only apparent, since it can be 
entirely determined by one more physical condition, namely 
the observer's reference frame. 

As a matter of fact the Weitzenbock and Riemann geome- 
tries are intrinsically related. Consider first the Cartan connec- 
tion ll24ll . r^v — e " /A e «v, defined in the Weitzenbock space- 
time. We can equally well write it as 

r^v = °r^ v + k /jXv > (2) 



where the T pXv are the Christoffel symbols, and is given 
by 

K /jXv = 2 (T^v + Tvl/j + TfiXv) ■ (3) 

is the contortion tensor defined in terms of the torsion 
tensor 7^ v , which is derived from r„x v . Indeed we have 

T\ v = d x e%-d v e a x . (4) 

The tetrad field transforms Lorentz indices into space-time 
ones and vice-versa, thus T pXv = e a/J T" Xv . 

The curvature of Cartan space-time vanishes identically, as 
it follows from substituting rU v in Eq. (fJJ, which yields the 
following identity 

eR(e) = -e{ l -T abc T abc + l -T' b <T bac - T a T a ) + 2d,{eT") , 

(5) 

where T p = T b and R(e) is the Riemannian scalar curva- 
ture constructed out of the tetrad field. 

The density lagrangian of General Relativity (GR), which 
leaves the field equations invariant under coordinate transfor- 
mations, is precisely the left hand-side of Eq. (0. Given that 
any divergence term in a lagrangian density makes no con- 
tribution to the field equations, we drop out the last term on 
the right-hand side of Eq. (|5]l and define the Teleparallel La- 
grangian density as 

L{e m ) = -Ke{^T" bc T ahc + -T abc T bac - T a T a ) - L M 

= -KeT." bc T ahc -L M , (6) 

where K—l/ (16ft), Lm is the Lagrangian density of matter 
fields, and T. abc is defined by the equality 

Yftbc ^abc -(- Y bac Y cab ^ -f- — (y\ ac T b T\ ab T c ^) (7) 

with T a = e a fi T> 1 . 

The equivalence between TEGR and GR is now easily un- 
derstood. Since the they have the same lagrangian density, the 
two theories yield the same dynamics. This does not mean that 
they predictions are identical. For instance, there is no quan- 
tity equivalent to the gravitational energy-momentum vector, 
which will be defined below, in the realm of General Relativ- 
ity. In fact every attempt to define the gravitational energy in 
the context of GR, from Einstein's equations, has failed. This 
approach has only lead to pseudo-tensors in GR. 

The field equations can be obtained from the variational dif- 
ferentiation of the Lagrangian density with respect to e afJ . Af- 
ter a few algebraic manipulations the resulting equality reads 

e aX e blJ d v (eZ bXv ) - e(T. bv a T bVfl - ~e afl T bcd Z bcd ) = ^eT afl , 

(8) 

where T clfJ = e a ^T^x = \ is the energy-momentum tensor 
of the matter fields. 

Explicit calculations show that Eq. (0 is equivalent to the 
Einstein equations j25tl . 
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The field equations can be rewritten as 

d v {el." lv ) = -;Ue%(^ + , (9) 
where v 41 is defined by the relation 

= K (4£^r 6c " - /"E^r M ) . (io) 

£«A,v j s a skew-symmetric tensor in the last two indices. 
This symmetry leads to the following local conservation law 



d x {et ak + eT aX ) = 



(11) 



from which it is possible to write down the continuity equation 

4- f d 3 xee a Jt°^ + T^) = - idS, \ee a + . 
at Jv Js 

We therefore interpret f^ A ' as the energy-momentum ten- 
sor of the gravitational field ll26ll . With this concept in mind 
we can define the total energy-momentum vector in a three- 
dimensional spatial volume V in a familiar way, as 



d xt 



(12) 



The above-defined energy-momentum vector is frame de- 
pendent and independent of the choice of coordinates. In or- 
der to analyze the features of gravitational field it is therefore 
mandatory to set up the reference frame first. There are count- 
less tetrad fields that match a given metric tensor. In other 
words, the same physical system can be seen under the optics 
of as many observers as desired. 

In the Hamiltonian formulation of TEGR the constraint 
equations are interpreted as energy, momentum and angular 
momentum equations for the gravitational field H27I1 . The 4- 
angular momentum of the gravitational field can be shown to 
have the expression 



r ab 



-4k / d s xe(Y: 
Jv 



aQb 



- T. b0a ) 



(13) 



The P a and L ab obey the algebra of the Poincar group II28II . 
The above definition, like Eq. (fT2l . is coordinate independent 
and changes with the observer. 

Although obtained in a Hamiltonian formulation, Eq. (fT3T i 
is readily interpre table. The angular momentum is, of course, 
the vector product between the momentum and coordinate. 
Since the tetrad fields are the dynamical variables, it is pos- 
sible to understand the meaning of Eq. ( fTTT i from the perspec- 
tive of the lagrangian formalism. The tetrad fields play the 
role of coordinates, when contracted with the total energy- 
momentum tensor, yields the above equation. 

As explained before, the choice of the tetrad field is not 
random. It is in fact intimatelly linked to the frame observer. 
For a given metric tensor, an infinity of possible frames exists, 
each of which completely characterized by the tetrad field. 
In order to fix the kinematical state of the observer in the 
three-dimensional space, we have to specify six components 
of the tetrad field, the other ten components being related to 



the metric tensor. To this end, we consider the acceleration 
tensor El H 



1 



>{0)ab + I a(0)b~ 1 b(0)a\ ! 



(14) 



here written in terms of the torsion tensor. 

Given a set of tetrad fields, the translational acceleration of 



the frame along a world-line C then follows from (j) 



(0) 



.(0 



and 



the angular velocity, from <|>(i)(n- Consequently, the acceler- 
ation tensor is a suitable candidate to geometrically describe 
an observer in space-time. It does not contain any dynamical 
features dependent on field equations and has been tested for 
teleparallel gravity in many situations I131T - 433I1 . 



III. THE SPACE-TIME OF A ROTATING NEUTRON STAR 

The most general form of the metric tensor describing the 
space-time generated by a configuration with axial symmetry 
is represented by the line element li34ll 

ds 2 = goodt 2 + g n dr 2 + g 22 dQ 2 + g^dty 2 + 2go3d<S?dt , (15) 

where all metric components depend on r and 8: g^ v — 
gfj V (r,Q). If g denotes the determinant of the metric tensor, 

we have that ^g = [gng22(g03g03 - goog3i)} 1/2 - 

For a rotating neutron star, with arbitrary angular velocity 
about the z-axis, we have the following components of the 
metric JUS: 

goo = - exp (y + p ) + r 2 co 2 sin 2 exp (y - p ) , 
g03 = -C0r 2 sin 2 9exp(y-p), 
gu =exp2oc, 
g22 = r 2 exp2a, 

g i3 = r 2 sin 2 9exp(y-p). (16) 



Here a, y, p and co are metric potentials, which are functions 
of r and 8. The matter inside the neutron star being described 
by a perfect fluid, thus the energy-momentum tensor has the 
form 



T" v = (e + p)U> J U v +pg< 



(17) 



where £, p and V are the energy density, the pressure and the 
fluid four-velocity field, respectively 113511 . Our assumption of 
axial symmetry makes V proportional to the time and angu- 
lar Killing vectors. Therefore, If (l,0,0,i2), where Q. is 
the angular velocity of the fluid, measured at infinity. In this 
model, to define the metric components in Eq. ( [ToT l, we have 
to specify an Equation of State (EOS) i37ll relating the energy 
density to the pressure and the angular velocity of the fluid. 

In order to analyze an axi-symmetrical spacetime within 
TEGR, we choose a stationary observer, who has to satisfy 
§ a b = 0. Adapted to this frame, the tetrad field reads 



( -A 

Jg\\ sin0cos(|) 

<Jgu sin0sin(|) 

w^iY cos 6 



-C \ 

/g22cos0cos(|) — Drsin9sin(|) 

/822 cos 6 sin (|) _Drsin6cos(|) 

-\/g22 sin 6 j 

(18) 
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The functions A,C and D are given by the equalities 



A = 
C = 

D = 



-goo) '/ 2 , 

g03 

(-goo) 1 / 2 



(r 2 sin 2 9)goo 



1/2 



(19) 



while 8 is defined by the relation 8 = go3g03 — goog33- 

We use a numerical method to solve the Einstein equa- 
tions producing the metric tensor. Our purpose is to es- 
tablish a relation between the features of the gravitational 
field, such as its angular momentum, and some experimen- 
tally measurable intrinsic attribute of the neutron star, such 
as the moment of inertia. Over the years, a convenient nu- 
merical method has been developed ll38ll . improved ||39"1 l40ll 
and modified fiUBll . We use the RNS code, available 
at |http: / /www. gravity .phys .uwm.edu/rns/| here modi- 
fied to deal with such teleparallel quantities as the energy- 
momentum vector and gravitational angular momentum. Our 
calculations use the non-dimensional quantities listed in 
Eqs. (4)-(13) of Ref. 39 and references therein. Here, how- 
ever we use \f~K, with K — c 4 /Ge c -, as the fundamental length 
scale of the system, instead of K' N I 2 , where K' is the poly- 
tropic constant, and /V is related to the adiabatic index, which 
would be more suitable to deal with poly tropic stars,. This as- 
sumes the matter to have no meridional motion, and the angu- 
lar velocity £2, as seen by an observer at rest at spatial infinity, 
to be constant. 

The program uses compact coordinates, /j and s, defined by 
the equality 



/j = cos 8 , 



R,. 



(20) 



where R e is r at the equator of the star. Thus s = 0.5 repre- 
sents a point on the equatorial surface of the star, while 5=1 
represents spatial infinity. 

Figure Q~]depicts an illustrative example of the output of the 
code. We have chosen the central energy e c /c 2 = 10 15 g/cm 3 
and adopted the equation of state in Ref. 46, here denoted by 
the acronym EOSA (dense neutron matter). The metric com- 
ponents indicated in the figure recover the Minkowski space- 
time in spherical coordinates at spatial infinity. 

The angular momentum and energy-momentum of the 
gravitational field will be discussed in the following sections. 



A. The Gravitational Angular Momentum of Neutron Stars 

This section uses the results developed in Ref. H3. Equa- 
tion (40) of that paper, which relates the z-component of the 
gravitational angular momentum to the metric, reads 



Kl)(2) _ 



-2k I 



\J-gao 



(21) 



Equation (f2Tb describes an axi-symmetrical spacetime, 
only the z-component of the angular momentum being 




FIG. 1 . Example of the metric behavior for fj — 0, using EOSA. Here 
Goo = goo, Gn = gli. G22 = g22, G33 = g33 and G03 = g03 • 10~ 3 - 



nonzero, and the analysis was restricted to slowly, rigidly ro- 
tating neutron stars. Under these conditions, the angular mo- 
mentum of the field was found to be proportional to the (very 
small) angular momentum of the source. Given that these con- 
ditions are unrealistic, we now want to study rapidly rotating 
stars in more complex configurations with a view to deter- 
mining how the field angular momentum depends on source 
features. We will use the shorthand L for LS l ^ 2 >, which will 
be dimensionless, unless otherwise stated. To transform back 
to the MKS system, all we have to do is to let L —> ^-^-L, 
where K is the fundamental length scale used in the numerical 
simulation. 

We now have a computational tool to simulate rapidly rotat- 
ing neutron stars and are able to investigate the behavior of the 
gravitational angular momentum. We will consider conditions 
that are more general than those in Ref. 53 

Figure [2] shows L as a function of s. The ratio between 
the polar and equatorial radii is 0.6. The curves peak approxi- 
mately at the surface of the star, then decay smoothly as s — > 1, 
a behavior in agreement with physical intuition, since we ex- 
pect L to be maximum near the surface of the star. 

The ratio between the polar radius (R p ) and the equatorial 
radius (R e ) defines the shape of the star. The larger this ratio, 
the smaller the gravitational angular momentum should be, 
since a unitary ratio, the highest possible value, corresponds 
to a spherical, static star. Figure [3] shows L as a function of 
the ratio R p /R e . The plots for EOSB, EOSC and EOSFPS 
are nearly coincident. In the intermediate regime, with R p /R e 
neither too small nor too close to unity, L is nearly linearly 
related to R p /R e . 

We are also interested in comparing the dependences of the 
field (L) and source (/) on the ratio R p /R e . We hence plot 
L/J as a function of R p /R e in Fig. [4] Again the curves for 
EOSB, EOSC and EOSFPS are superimposed, a feature also 
found in Figs. [5] and [6] Figure|4]shows that in rapidly rotating 
neutron stars the gravitational angular momentum is compara- 
ble to the angular momentum of the source. The gravitational 




FIG. 2. Spatial distributions of the gravitational angular mo- 
mentum for the indicated equations of state. The central energy 
is E c /c 2 = 10 I5 g/cm 3 and the equations of state are as follows: 
EOSA [46] (dense neutron matter), EOSB [48] (hyperonic mat- 
ter), EOSBBB1 1 49] (asymmetric nuclear matter, derived from the 
Brueckner-Bethe-Goldstone many-body theory), EOSC [50] (dense 
hvnernnic matter) and F.O.SFP.S 1511 f neutron matter for an imnrnved 
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FIG. 3. Gravitational angular momentum as a function of the ratio 
between the polar and equatorial radii. The legend follows the con- 
vention of Fig. [2] 
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FIG. 4. Ratio between the gravitational angular momentum and the 
angular momentum of the source at the spatial infinity as a function 
of the ratio between the polar and equatorial radii. The legend fol- 
lows the convention of Fig. [2] 



ure [5] shows that the gravitational angular momentum grows 
linearly with the moment of inertia, except at the edges of 
the figure. This is what we expect on intuitive grounds. The 




field can therefore be detected, since the angular momentum 
of the source is experimentally accessible. We have found 
that the ratio L/J is proportional to the inverse of R p /R e . For 
a given stable star, L is therefore proportional to R e /R p J. 

Our computations make no reference to binary systems. 
Nonetheless, the conclusion that L and J have comparable 
magnitudes shows that the apparent non-conservation of an- 
gular momentum in such systems is not necessarily linked to 
the intensity of the emitted gravitational waves. While beyond 
the scope of this paper, additional investigation along this line 
of reasoning seems undoubtedly interesting. 

Figures (0 and d6) present numerical results for the depen- 
dence of L on two other parameters of the star, the moment of 
inertia and the angular velocity of the fluid, respectively. Fig- 



FIG. 5. Gravitational angular momentum as a function of the mo- 
ment of inertia. The legend follows the convention of Fig. [2] 

In Fig. [6] by contrast, the relation between the gravitational 
angular momentum and the angular velocity of the fluid is 
non-linear. This is also reasonable, since the angular momen- 
tum varies as the angular velocity of the source grows, because 
the rotation is not rigid. 

B. The Energy-Momentum Vector of Neutron Stars 

For completeness, we now present numerical results for the 
energy and momentum. Given that the metric is stationary, 
it is reasonable to expect no energy flux. Even if there is an 
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FIG. 6. Gravitational angular momentum as a function of the angular 
velocity of the fluid. The legend follows the convention of Fig. [2] 



initial flux, the star will emit gravitational waves until stability 
is reached, at which point our method becomes trustworthy. 
The combination of the tetrad field in Eq. ([T8l > with Eq. (0 
leads, after some algebraic manipulation, to the result 



4eE (0)01 = 



1 



2^11^228-8 



\/(-goo8) giign 
■ V 8228 1\ - {822) 3/2 Vgiigoi ( - )r 
+ 2y/-go Q (g22) 3/2 V8giisin9 + 
+ 800 



V dr 

d80i 



n3/2 , f^lA 

122) 1 \/eu I -^r— I 



V dr 



(22) 



Next, on the basis of Eq. ( fT2b . we examine the energy as a 
function of s. We work with dimensionless quantities and let 
E denote P^O). To recover the energy from the dimensionless 
variable, we let E -> (c 4 y/K/G)E, where y/K is the funda- 
mental length scale. 



First, a simple check on the code. Figure Q shows the total energy, i. e., the energy in a hyper-surface of infinite radius, as a 
function of the mass. Although straight over a large fraction of the displayed range, the solid lines in the first five panels bend up- 
or downwards at the highest masses. That the curvatures signal the breakdown of the nonrelativistic equations of state is shown 
by the rectilinear plot in the last panel of Fig. [7] which contrasts with the other panels because the EOS in Ref. ll52ll includes 
relativistic corrections. 
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FIG. 7. Gravitational energy as a function of the total mass at the spatial infinity. In this figure m = 10 M, where M is the total mass of the 
star. The last panel, labeled EOSAPR, has been constructed from the equation of state in Ref. 1 52], whci includes relativistic corrections. 



Figure <(8j shows the ratio between the energy and mass as 
a function of s. For each equation of state, E/M rises up to 
the surface of the star and then decreases towards a constant 
as s — > 1, showing that the energy at infinity is proportional to 
the mass — a conclusion that cannot be directly inferred from 
Fig. [7] The slow decay beyond s rj 0.45 shows that the rota- 
tional energy is very small. 

Fig. [9] shows the gravitational energy as a function of the 
squared angular velocity of the source for each of the indi- 



cated equations of state. As the angular velocity grows, the 
initially straight plots bend downards. While, as already ex- 
plained, the deviations from linearity in Fig.|7]reflect the inac- 
curacy of the non-relativistic equations of state, those in Fig. [9] 
have physical origin: they indicate that the moment of inertia, 
which relates the energy of the field to the angular velocity, 
depends on the angular velocity of the fluid. 

We have also computed the momentum of the gravitational 
field, using the spatial indices in Eq. (TTZl ). The relevant com- 
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FIG. 8. Spatial distribution of the ratio between the gravitational 
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and E^ 01 = 0. 

The momentum is therefore zero, a result that can only 
be mathematically understood if we recall that the right-hand 
side of Eq. (fTSl i has been integrated over the coordinate (j). 
Physically, we expect the rotating star to transfer no momen- 
tum to the field. 

IV. CONCLUSION 



We have analyzed the physics of neutron stars from the point 
of view of teleparallel gravity. We have determined the energy 
and angular momentum as functions of the compact spatial 
coordinate. In particular, we have computed the ratio between 
the field and source angular momenta. For rapidly rotating 
neutron stars the magnitudes of the two angular momenta are 
comparable. Figure [6] showed that, as expected, the ratio van- 
ishes in the limit of slow rotations. As shown by Fig. [6] at 
high angular velocities the field angular momentum ceases to 
be proportional to the source angular momentum, an indica- 
tion that the moment of inertia has become strongly dependent 
on the angular velocity of the field. 

We have calculated the ratio between the gravitational en- 
ergy and the mass as a function of the compact coordinate 
s. This ratio is maximized near the surface of the star. Be- 
yond the surface, it decays to a constant at infinity, a value 
that is weakly dependent on the equation of state. We are 
lead to conclude that, as one would expect, each type of neu- 
tron star is characterized by a universal linear relation between 
the gravitational energy and mass at spatial infinity. Although 
we have found small deviations from linearity for sufficiently 
large masses, the comparison in Fig. |7]has linked those devia- 
tions to inaccuracies in the non-relativistic equations of state. 



FIG. 9. Gravitational energy as a function of the squared angular 
velocity of the source. The legend follows the convention of Fig. [2] 



ponents are 
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